Frequency and damping of hydrodynamic modes in a trapped Bose-condensed gas 
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Recently it was shown that the Landau-Khalatniltov two-fluid hydrodynamics describes the 
collision-dominated region of a trapped Bose condensate interacting with a thermal cloud. We 
use these equations to discuss the low frequency hydrodynamic collective modes in a trapped Bose 
gas at finite temperatures. We derive a variational expressions based on these equations for both 
the frequency and damping of collective modes. A new feature is our use of frequency-dependent 
transport coefficients, which produce a natural cutoff by eliminating the coUisionless low-density tail 
of the thermal cloud. Above the superfluid transition, our expression for the damping in trapped 
inhomogeneous gases is analogous to the result first obtained by Landau and Lifshitz for uniform 
classical fluids. We also use the moment method to discuss the crossover from the coUisionless to the 
hydrodynamic region. Recent data for the monopole-quadrupole mode in the hydrodynamic region 
of a trapped gas of metastable "^He is discussed. We also present calculations for the damping of 
the analogous m — monopole-quadrupole condensate mode in the superfluid phase. 



I. INTRODUCTION 



^ • PACS numbers: 03.75.Kk, 05.30.Jp 

B 

. 

^ , In recent papers, Zaremba and the authors have derived a closed set of the two-fluid hydrodynamic equation of 
Q ' a trapped Bose-condensed gas starting from a simplified microscopic model describing the coupled dynamics of the 
O condensate and noncondensate atoms 0, Qi S S 11 These equations can be written in the Landau-Khalatnikov (LK) 
. form, well known in the study of superfiuid ''He |a, Q . These simplified hydrodynamic equations include dissipative 
I ' terms associated with the shear viscosity, the thermal conductivity, and the four second- viscosity coefficients. Explicit 
^ , formulas for these transport coefficients were obtained in Ref. ^ and used to define three characteristic transport 
^\ relaxation times These define the crossover between the coUisionless and hydrodynamic regions. Detailed calcu- 
I lations of these transport relaxation times in a trapped Bose gas Q shows that the collisions between the condensate 
• and noncondensate enhances the transport relaxation rates significantly in the MIT data @ , so that one is in the hy- 
' drodynamic region below Tbec- We also note that the recent Bose condensate observed in metastable He* [l0lllTlll2 | 

■ appears to be well within the collision-dominated hydrodynamic region, even above the Bose-Einstein condensation 
' temperature Tbec- This is because of the relatively large density of He* atoms and their large s-wave scattering 
', length. 

■ In the present paper, we derive a general expression for the frequency and damping of hydrodynamic collective 
modes in a trapped Bose-condensed gas at finite temperatures, starting from the two-fiuid hydrodynamic equations 
derived in Ref. ;5]. These two-fluid equations are briefly reviewed in Sectional and reformulated as a closed set of 

' O [ equations for the condensate and noncondensate velocity flelds. In Section IlIII we derive a variational expression for 
^ ■ undamped normal-mode frequencies in the Landau limit {ujt 1), extending an approach flrst developed in Ref. 0]. 
In Section Hvl we obtain a general expression for the damping, which only depends on knowing the undamped normal- 
mode solutions. This kind of expression is very convenient in working out the damping of hydrodynamic modes in 
trapped Bose gases, as first pointed out by Kavoulakis et al. |0. As an illustration, in Section IVll we give a detailed 
discussion of the the m — monopole-quadrupole collective mode above Tbec studied in the recent experiments [l2l| . 
In section rVlII we also calculate the damping of the m = hydrodynamic mode in the superfluid phase. 
^ Appendix A gives some details of the damping calculations based on the use of frequency-dependent transport 
coefficients. The moment method for a degenerate normal Bose gas is reviewed in Appendix B. 

II. TWO-FLUID HYDRODYNAMICS OF A TRAPPED BOSE GAS: A REVIEW 

In this paper, we consider a Bose-condensed gas confined in an external anisotropic harmonic trap potential 

t/ext(r) = |(..Ix2+c.y-|-c.,V), (1) 
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Our starting point is the two-fluid hydrodynamics derived from the finite-temperature kinetic theory by Zaremba, 
Nikuni, and Griffin (ZNG) j3j|]. In the ZNG theory, the coupled dynamics of the condensate and noncondensate Jjj is 
described by the generahzed Gross-Pitaevskii (GP) equation for the condensate wavefunction <I>(r,t) 



ih 



a$(r, t) 
dt 



2m 



+ C/cxt(r) -t- gndr, t) + 2.gn(r, t) - ii?(r, t) 



and the semi-classical kinetic equation for the noncondensate distribution function /(r,p,t) 
3/(r,p,i) , p 



dt 



Vrf{r,p,t) ~ VU{r,t) ■ Vp/(r,p,t) = C,2[f, $] + C22[/]. 



Here nc(r,t) = |$(r,t)p is the condensate density, and n{r,t) is the noncondensate density. 



n{r, t) 



dp 



(2nhy 



:/(r,P,i), 



(2) 



(3) 



(4) 



and U{r,t) — t/ext(r) -I- 2g[nc{r,t) + n{r,t)] is the time-dependent effective potential acting on the noncondensate, 
including the Hartree-Fock (HF) mean field. As usual, we approximate the interaction in the s-wave scattering 
approximation g = Anfi^a/m. The dissipative term i?(r,t) in the generalized GP equation (0) is due to the coUisional 
exchange of atoms in the condensate and noncondensate. This is related to the C12 collision integral in iPJ, namely 



Rir,t) 



ri2(r,i)= I ^^^Ci2[/(r,p,i),$(r,t)]. 



(5) 



Explicit expressions for the two collision integrals (C22 and C12) in the kinetic equation ^ can be found in Rcf. 3]. 

The GP equation (5) can be written in the hydrodynamic form in terms of the amplitude and phase of $(r,i) = 
Vnc(r, i)e**''''*), which leads to 



duc 

at 

d_ 

dt 



V-(n,v,) = -ri2[/,$] 



(6a) 
(6b) 



where the superfluid velocity is Vc = h^6(r, t)/m and the condensate chemical potential is given by 



Mc(r,t) 



2m^Jnc{r,t) 



+ C/cxt(r) -I- gnc(r, t) -f 2gh{Y, t) . 



(7) 



One sees that ri2 in Eq. (|6a|l plays the role of a "source function" in the continuity equation for the condensate, 
arising from the fact that C12 collisions do not conserve the number of condensate atoms 3]. 

Hydrodynamic equations for the noncondensate can be derived by following the standard procedure first developed 
in the kinetic theory of classical gases. We take moments of the kinetic equation Q with respect to l,p and to 
derive the most general form of "hydrodynamic-type equations" for the non-condensate. These moment equations 
take the form (/i and v are Cartesian components): 



^ + V(nv„)=ri2[/], 



dt 



'-'nfj. 



dP^, .dU 



de 

— -h V • (ev„) = - V • Q - D^^P^^ 



ri2[/]. 



(8a) 
(8b) 

(8c) 



Here and elsewhere, repeated Greek subscripts are summed. The noncondensate density was defined earlier in Eq. Q), 
while the noncondensate local velocity v„(r, t) is defined by 



n(r,i)v„(r, t) 



dp P 

{2Trh)^ m 



/(r,P,t). 



(9) 
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In addition, we have introduced the following quantities: 



P,Ar,t) ^ mj ^(^^.„,) (^_^/(r,p,i), (10a) 



Finally, the symmetric rate-of-strain tensor appearing in Eq. (|8c|l is defined as 

D^^^{^,t) = 7^\^ + ^] ■ (11) 



The "hydrodynamic" equations are exact, but not closed as they stand. 

We next apply the Chapman-Enskog procedure to obtain a closed set of hydrodynamic equations. This procedure 
yields ^ 



= ^M^-P - 2?y l^D,,, - ^TrDS^^j , (12) 

Q = -kVT, (13) 

where P is the kinetic pressure and T is the temperature. The above formulas involve the position-dependent shear 
viscosity rj and the thermal conductivity k. These local position-dependent transport coefficients will be discussed in 
more detail below. 

In order to study small amplitude oscillations, we linearize the hydrodynamic equations around static equilibrium 
as Uc = rico + Sric, Vc = (5vc, n = hq + 5h, v„ = (5v„, P ^ Pq + 6P, where the subscript denotes static equilibrium. 
The equilibrium condensate density profile is determine by (within the Thomas-Fermi approximation, which neglects 
the quantum pressure or kinetic energy term in the GP equation) 

nco(r) = -[^ico - C/oxt(r)] - 2no(r), (14) 
9 

while the equilibrium noncondcnsate distribution function /o is given by 

/o(r,p) = -— T— 1^ ■ . (15) 

The local density fiQ (r) and the local kinetic pressure Pq (r) of the noncondensate atoms are given from Eq. (fTK|l as 

^o(r) ^.g3/2(^o), (16) 

^o('-) = ^55/2(^0), (17) 

where ^o(r) = e^'^''°~^°'^^^^/''^'^'' is the local equilibrium fugacity, A = {2'n:h^ /mk^T^))^/'^ is the thermal de Broglie 
wavelength, and gn{z) = X^z^i -^V"- 

The linearized hydrodynamic equations for the condensate are given by 

^ = - V • (n,o<5v,) - 5ri2, (18a) 
at 

m-^ = -gV5{5n^ + 25n), (18b) 
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while the Hnearized hydrodynamic equations for the noncondensate atoms are given by 



dSn 



mnQ 



V • (no(5v„) = STi2 
dSP 



dt 

dSP 5 



f.~dUo ^ . dSn d . 

"a """a 25^0-^ h — <^ 2?7 

dx^, dx^ dXf, dx, 



dt 



Uo)STi2 + • (kVST). 



(19a) 
(19b) 

(19c) 



The above equations involve the fluctuations of the source function (5ri2 and the temperature ST. These can also be 
written in terms of the condensate and noncondensate velocity fields Q. One finds ^ri2 = ^T^^ + ^r^2\ where 



(1) 

2 



6T 



(2) 



V 

d 



12 

dST 

~dr 



(1) Scthcti 



1 



V • (kVST), 



T^O'l'ji 12 ■ 



(20) 
(21) 
(22) 



Here an and cti only involve the static local thermodynamic functions of the gas, and are defined in Eqs. (25) and 
(57) of Ref. [|. In Eq. r,, is a new relaxation time describing how fast the condensate and noncondensate atoms 
reach diffusive local equilibrium with each other. More explicitly, it is given by 



(23) 



where ri2 is a collision time of the condensate atoms with the noncondensate atoms as defined in Ref. [Bl • 

The above two-fluid hydrodynamic equations have dissipative terms involving the shear viscosity rj and the thermal 
conductivity k. These transport coefficients are given by the following expressions 



m I 2g3/2(zo) 



55/2(20) 



53/2(20) 



(24) 



77 = r^nofceTo 



55/2(20) 
.93/2(20) 



(25) 



where zg is the local fugacity defined earlier. These expressions involve the characteristic transport relaxation times 
and T,j, which are defined in Rcfs. 5, "s]. The three relaxation times t^^jT^ and characterize how fast the 
two-component system reaches local equilibrium, and thus they define the crossover frequency between the collision- 
dominated hydrodynamic region and the so-called coUisionless region dominated by mean fields. We note that the 
two- fluid equations 118(1 and 1)19(1 can also be rewritten in the LK form 0| , which involve the four second viscosity 
coefficients C1JC27C3J and ^4- As shown in Ref. ^5], they are all related to the relaxation time Tfj, as follows: 



Ci 



C4 



3m 



. 5^ . 5 

C2 = -^Cr^T^, C3 = 



(26) 



To derive a closed set of equations for the velocity fields Vc and v„, we take time derivatives of l|6b(l and (|19b|l . We 

obtain 



9^ 



= .gV[V • (neo<5vc) + 2V • (no'5v„)] - gV5Yi2. 



(27) 



and 



5 



d 

+2gnQ^-—[V ■ {nco5^fc) + V • (no(5v„)] 

O KJJbn O UJUn 



1 



dD 



We then look for normal-mode solutions of (|27|l and (|28(l of the form 

v„(r, t) = u„(r)e-*"*, v,(r, i) = u,(r)e-*"*. 
In this case, the coupled equations (|27|l and (|28l) reduce to 

mw^Uc = -5V[V • (ricoUc)] - 2gV[V • (nou„)] + 5V5ri2,i^[u„, u^], 



(28) 

(29) 
(30) 



mriQUJ Unn 



-It^IV • (Pou„)] + ^^(u„ . VPo) - V . (nou„)|^ 



3 dxi. 
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d 



2ghQ^—[V ■ (ncoUc) + V • (nou„)] 

IdUo 2 9(5ri2,c^[u„,Uc] 

-dri2,i^[u„,UcJ - -guco- 



3 



fen 



The symmetric tensor u^^ is defined by 



iiu ~ {y ■ u„) 



1 / dUniy . aMj, 



2 d 



2 y dxn dxu 



(31) 



(32) 



The source function 5Ti2 appearing in the above equations can be expressed in terms of of u„ and Uc as 5Ti2 = 
<5ri2,c^[u„,Uc]e""^*, where 



^ri2,i^[u„, Uc] = (5r^y [U„, Uc] + (5ri2^^[u„, Uc], 



(2) 



with 



Sv'^y} [u„, Uc] = o-ff <{ V • [nco(uc - u„)] + |ncoV • u„ 



(33) 
(34) 



^^f2u> = «WTp<5r^2^[u„,Uc] - ^^4^V • (kJT^[u„,Uc 

Similarly, the temperature fluctuation is given by 5T — iST^jlun-, Uc], where 



(5r„[u„,Uc] = - 



a; 



4to(V • u„) + ^ai,5ri'2)[u„,uc 
3 3no 



(35) 



(36) 



Using these results in Eqs. ipUjl and (PJ), we see that we have obtained a closed set of equations for both local velocity 
components u„ and Uc. 



III. UNDAMPED NORMAL MODE FREQUENCY 



We first consider the undamped normal-mode solutions of our hydrodynamic equations, neglecting all hydrodynamic 
dissipation. Formally this means that we take 77, k, ^ in the two-fluid hydrodynamic equations. As discussed 
in Refs. i5i], this limit corresponds to the Landau two-fluid hydrodynamics without dissipation. In this limit, the 
coupled equations for u„ and Uc simplify to 

mio'^Uc = -.gV[V-(ncoUc)]-2.gV[V-(noU„)]+5V(5rj2^^^[u„,Uc], (37) 

muj^Un = - 3 V[V • (PoU„)] + 3 V[u„ • VPo] - 3 V (^nco^r^y^^iu,,, Uc]) 
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-5r^2^[u„, Uc]VUo - V • (noU„)Vf/o - 2.gnoV[V • (noU„) + V • (^coUc)], (38) 

In general, solutions of these coupled hydrodynamic equations are very complicated. However, one can reformulate 
this problem so that the solutions are given in terms of a variational functional. The present discussion closely follows 
the variational analysis developed in Ref. ^ for two-fluid hydrodynamic equations, introduced in Ref. Jj, which 
omitted the contribution of the source term from C12 collisions, i.e., ST12 = 0. However, as Ref. Q showed, i5ri2 
plays a crucial role in obtaining the correct Landau two-fluid hydrodynamic limit. The same formalism was also used 
in Ref. Q to calculated the hydrodynamic mode frequencies for a trapped Bose gas above Tbec- 
Introducing the six-component local velocity vector 



(39) 



we combine the coupled equations for u„ and Uc in (|38() and H37() into a matrix equation 

Lu = Lu^Du. (40) 



The 6x6 matrix L has the block structure, 



with the 3x3 matrix elements being defined as 



Lii L12 
L21 L22 



(41) 



(Liiu„)^ = _^_^[V.(Pou„)]-f ^^(u„- VPo) 

-V • (nou)- 2gno^ — [V • (nou„)] 

(|9?^cO<5r«[u„,0]) + ^rW[u„,0]^, (42a) 



{Li2VLc)i, = -2.gno^^V • (ricoUc) - 0gnco<^r^y [0, u^] 

+Jriy[0,ujg, ' (42b) 

(L2iu„)^ = -~2gnco-^'^ ■ {fioxin) + gncQ^^5T^y}[\in,0], (42c) 

{L22Uc)f, = -gricQ-^V ■ {ticqUc) + gncQ-^ST[^^[0,Uc]. (42d) 

Similarly, the matrix D is block-diagonal (I?i2 = £'21 = 0), with elements Dn — mfiQl, D22 — mricoi- Here, 1 is a 
3x3 unit matrix. We note that the matrix L has the Hermitian property 

J drW ■ (Lu) = y" firu • {Lu'). (43) 

The coupled equations (|37|l and H38|l can be rewritten in terms of the variational functional 

J[u„,Ue] ^ , (44) 

K [u„, UcJ 

where 

U[un, ^c] = \j dru ■ (Lu), K[un, ^c] = lj dru ■ (Du). (45) 

Using the Hermitian property of L in Eq. H43() , one can prove that the requirement that the functional J be stationary 
leads to the required equations in Eqs. (|30|l and H31(l and is identified with the stationary value of the functional J. 
One can therefore evaluate the collective mode frequency using a variational ansatz for u„ and in the variational 
functional J[u„,Uc]. 
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IV. GENERAL EXPRESSION FOR DAMPING OF HYDRODYNAMIC MODES 

In this section, we derive a general expression of hydrodynamic damping of a collective mode due to transport 
coefficients. The general expression for hydrodynamic damping of collective modes in a trapped Bose gas was first 
derived in Ref. [13 above Tbec- 

We now include hydrodynamic dissipation involving k, 77 and in the two-fluid equations. Similarly to Eq. H40|l . 
one can write the coupled equations for u„ and u^, which are given in Eqs. (|30|l and (|31|) . in a matrix form 

Lu + Fu = uj^Du. (46) 

Here F represents the dissipative terms in the two-fluid equations and has the block structure 



with the matrix elements 



Fu F12 
F21 F22 



ic^r^(5r«[u„,0] - ^^^V • (KV5r^[u„,0]) 



(47) 



(i^iiu„)^ = -— ^-gn,o{^uJT^)ST['J[n,,,0] - i^- + ^'^J^ ] V • (AcV<5r^[u„, 0]) 



'f^V.(.VJr.[u„,0])jf^ 

3gno J dXfj, 

d ^ f 1, 



-iuj^—2'q ( Uf,^ - -V • Un6^^ ] (48a) 



iLUT^SrY2'[0,u,] ~ • (kV^T^[0,uJ)| (48b) 



d 



d 



^c^T^<5^(y[u„,0] - ^J^^V • (^V5T^[u„,0]) 



(48c) 
(48d) 



In our subsequent analysis, we treat as a small perturbation to the undamped equation in Eq. (|40|) . To flnd a 
solution for u„,Uc including damping, we expand the vector u in terms of undamped normal- mode solutions (a 
is the mode index): 

U=^CaUa, LWa=^\DWa- (49) 

Q 

From the Hermitian property of the operator L, one can show that these normal-mode solutions satisfy the orthonor- 
mality relation ^ 

J drua- (Dufj) = Sap ■ (50) 

We note that in the above relation, we assume that the normal- mode solutions are normalized. Making use of ((SUJ, 
we obtain a linear equation for the coefficient C'a '■ 

(lu^ - ujl)Ca = J dvua ■ J2 Ca'FvLa' = J2 y^o^'Ca', (51) 

a' a' 

where the matrix element Vaa' is defined by 

Vaa' = J drua ■ Fu^' ■ (52) 
We note that Vaa' also depend on the frequency lu. 
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Expanding the mode frequency to first order in the perturbation F as w = uja + Acua , we find 

AWq = -^VaaU^a-o, (53) 



where 



We thus find that AcJq — — iFc, where the damping rate Fq of the mode a is given by 



V-u„„j }. (54) 



J \ O-H Jo 



-\ -1 

2 / drm(ncoMcQ + ^o^L) 



(55) 

Here we exphcitly display the normaUzation factor in Eq. (|55ll . This expression for the damping rate can also be 
written in terms of the second viscosity coefficients: 

Ta = y"dr|c2(V • u„q)^ + 2Ci(V • u„a)V • [mnco(ucQ - u„a)] 



+C3{V • [mUrMiUca - U„q)]}^ + 7^ | V(5T^„ [u„a , Uca] ^ 



To' 

2~ 



-7: 1 7^ -^^uv^ ' Una 

aXfj. OXv 3 

-1 



2/dr^(n.o.L + -oO 



(56) 



The formula in Eq. (|56|) for the damping rate can be understood in terms of the entrow production 16]. The local 
entropy production rate Rg (r , t) in the two-fiuid equations ^ is given in Eq. (87) of Ref. j5j . Assuming a normal- mode 
oscillation of the form 

V„(r,t) = Una{-r) COS. UJ at, Vc(r, i) = Uca(r)cosWai, (57) 
we find that the time average of the total entropy production rate is given by 

= ^ y dr|c2(V • u„a)2 + 2Ci(V • u„a)V • [mncoiuca - u„q)] 

+C3{V • [mricoiUca - U„„)]}2 + ^\yST^JUna,Uca]\'^ 

-to 

+H^ + ^~3^-^-H /■ 

On the other hand, the total mechanical energy is 

(£^moch) = 7: / dr(mrioM^„ + mncou^a). (59) 



(60) 



One can then write the damping rate Ta as 

2(-£'mcch) 
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This general expression for damping was first given in the classic work in Landau and Lifshitz (LL) for classical 
fluids. It was later used by Kavoulakis et al [l^l to study damping in trapped Bose gases above Tbec- This kind of 
LL damping formula was discussed in the case of superfluid "^He by Wilks |3| ■ 

So far we have not dealt with the problem arising from the fact that in a trapped Bose gas, the decreasing density 
in the tail of the thermal cloud always leads the breakdown of the hydrodynamic description. As pointed out by 
Kavoulakis et al. llJ. this causes trouble in using Eq. (|55|l or H56(l to evaluate the damping of modes in a trapped 
Bose gas. In Refs. |l3lll5j . this problem was handled in a physically motivated but ad hoc manner, by introducing a 
spatial cutoff in the integral. 

In this paper, we propose a new, more microscopic, procedure to deal with this problem. As we discussed in Ref. |^, 
the fact that the condensate and noncondensate atoms are not in complete local equilibrium can be taken into account 
by introducing the frequency-dependent second viscosity coefficients 

c.(^) = r^^- (61) 

Similarly, one can also introduce the frequency dependence in the shear viscosity and the thermal conductivity as 

= r;H = -^^. (62) 

1 — lUJT^ 1 — lUJTri 

In the Appendix A, we give a more detailed discussion and derivation of these frequency-dependent transport coef- 
ficients starting from the kinetic equation. Replacing the transport coefficients in Eq. I|56|) with k,{lj a) i a) and 
C,i{uja), and taking real part, we find the damping rate of a collective mode in a trapped Bose gas 

J \l + {LOaT^,Y OH ^ ' l + (w„T„)-^ro 



1 ri ( dUnau ^ dUnafj, 2^ u (5 ^ ^ 

f + (UaTr,)^ 2 \ dx^ dXy 3 



n -1 

2 I ~ „,2 



2 / drm(7icoMcQ + "ou„q) 



(63) 



We recall that (5T^„ [una, Mcq] and 5T''-^2[unaiUco\ are defined in Eas. H34l) and l|36|) in terms of the velocities Una and 
Uca- This result in Eq. (|63|) allows us to calculate the hydrodynamic damping due to various transport processes in 
a trapped Bose gas (both above and below Tbec) and it is is the major new result of this paper. The frequency- 
dependent transport coefficients in Eq. H63|) automatically yield the factors 1/[1 + (waT^)]^, (i = fi,K,ri), which 
effectively introduce a spatial cutoff for ujaTi > I, i.e., when hydrodynamics breaks down and we enter the coUisionless 
regime. 

V. A UNIFORM BOSE GAS 

As an illustration of the physics implied by our results in Section [IVI we study first and second sound in a uniform 
Bose-condensed gas using our variational expressions for the frequency and damping. In a dilute gas, first sound 
mainly involves the noncondensate oscillation (u„ ^ Uc), while second sound mainly involve the condensate oscillation 
(uc S> u„) (see, for example, Ref. To a first approximation, we can simply use Uc = for first sound and u„ = 

for second sound. Using the plane- wave solution u„, Uc oc fc cosk • r in the variational formulas, we find 

uji = Uik - iVi (i = 1, 2), (64) 

where the two sound velocities are given by 



^1 — Z — ~r — Z 

3mno ^ %mh{) 



(65) 



^""(1--h), (66) 



TO 
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and the damping rates are given by 

fc2 



2mno 



4 



mUcoiCi + C4) 



9ul 



3ho J 



(67) 



r, = — 



AkToHco 2 



(68) 



These results agree with Eqs. (C11)-(C14) of Appendix C in Ref. 0], which were directly derived from the Landau- 
Khalatnikov two-fluid equations for a dilute Bose gas. We note that the two sound velocities ui and U2 are slightly 
different from those given in Ref. 0|, because Ref. ^3 worked in the limit wr^ ^ 1, and neglected the source term 
5T 12- As shown in Ref. 0, however, these differences are quantitatively very small in the case of a weakly-interacting 
Bose gas. 



VI. MONOPOLE-QUADRUPOLE MODE IN A DEGENERATE NORMAL BOSE GAS 

Let us now consider collective modes in a trapped noncondensed Bose gas above Tbec- Here we consider the m = 
monopole-quadrupole collective mode in an an axi-symmetric trap {ux = Uy = ui± ^ Uz)- This type of collective 
mode above Tbec was first observed in the pioneering MIT experiment 0, but the density was not large enough to 
probe the hydrodynamic regime (see, however, discussions in Ref. Q). More rece ntly , however, ENS experiments with 
metastable He* atoms studied the m = mode in a high-density thermal cloud (laj . 

Above Tbec: the Hartree-Fock mean field is negligible and thus the equilibrium density is simply given by Eq. I|16|) 
with zo(r) = e'^I'^""^"'''^'"^!. The chemical potential /ig is determined as a function of the temperature through 
N = (k^T / TiujY' g^i^zo) ^ where Co = {ujx^y^zY^^ ■ The hydrodynamic modes in a trapped Bose gas with the to = 
symmetry were first discussed by Griffin, Wu, and Stringari |l8|| . The two normal- mode frequencies are temperature- 
independent, and are given by |l8| 

[5wi -I- Awl ± a/25u;4 +I6wf-32w2tj2] . (69) 

o • 

The corresponding velocity field is given by 

u„ = {ax,ay,bz), (70) 
where the coefficients a and b satisfy the following relations: 

Kavoulakis et al. discussed the hydrodynamic damping of this to = mode using the LL formula, with a 

spatial cutoff to deal with the crossover from the hydrodynamic to coUisionless regime in the tail of the thermal cloud. 
In contrast, we calculate the damping rate using Eq. (|63|l . which eliminate the coUisionless regime in the tail of the 
thermal cloud through the use of the frequency-dependent transport coefficients. Above Tbec, there is no contribution 
from the second viscosity transport coefficients. Moreover, in the to = mode above TbeCj one can show VT = 
[TsLITsI l and thus the thermal conductivity makes no contribution to the damping of the to = mode. Thus, only the 
shear viscosity in Eq. (|63l) contributes to the damping of this low- frequency collective mode. Using Eq. ((TDJ, we find 
(see also Ref. jU]) 



j dr [a'i_{x^ + y^) + h'^_z'^]mha{Y) 



Making use of the fact that the equilibrium density profile given by Eq. H16|) is a function of U^xt (r) , one can rewrite 
Eq. lO as 



3' J l + [f^-r,(r)]2 ^^^^ 
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In Eq. (|73() . the factor involving the coefRcients a_ and 6_ can be written in a simple form in terms of the undamped 
frequencies r2_ and f2+ as 



(74) 



Here we have used Eq. (|71f) and Eq. (|69|) . We thus obtain the following simple expression for the damping rate 



2(n2__n2_) 

where we have introduced a new relaxation time [see Eq. (|25|) ]: 

^o(r)r,(r) 



(fii -4a;2)(r!l -4a;l), (75) 



dr- 



TO 

"3 



Here we have used the relation [see Eqs. (|16ll . (|17|) and (|25() ] r/(r) = _Po(r)'n,, (r). Using the equilibrium relation 
VPq + "-oVJ/oxt = 0, which is valid when neglecting the HE mean field, one can reduce Eq. (|76|l to 



dr 



1 + {n^T, 



)2 

-• (77) 

rfrPo(r) 

Before presentin g re sults given by an explicit evaluation of the relaxation time f in Eq. (|76() using the parameters 
for the ENS trap jl^l, it is useful to comment on the relation between the present calculation of hydrodynamic 
damping and the moment method developed by Guery-Odelin et al. . These authors applied the moment method 
to the classical-gas Boltzmann equation. It is straightforward to generalize their method to the kinetic equation for 
a Bose-degenerate gas (see, for example, Ref. |20j|l. This moment method is briefly reviewed in Appendix B, and 
here we simply give to the final results. In the moment method, collisions are characterized by a single parameter, 
the quadrupole relaxation time r defined in Eq. I|B7|) . In the hydrodynamic limit cj^t <C 1, the moment method 
reproduces the hydrodynamic frequency in Eq. 169() , but the damping is now given by Eq. ljB12|l . The moment result 
for the damping has the same form as Eq. (|75|) . except that r replaces f. Both f and r arc related to the same 
position-dependent viscous relaxation time Tv,(r), but involve different spatial averages [see Eq. (I76|) and Eq. (|B7p ]. 

Our evaluation of the damping F in Eq. (|75|l is based on the ENS trap parameters 0, ll^J: trap frequencies 
LLi\^l2'n — 988,ct;2/27r 115, total number of atoms N = 8.2 x 10^, and s-wave scattering length a — 16nm. The 
ideal Bose gas transition temperature is given by Tc = (?i[I'/fcB)(Af/1.202)^/'^ = 4.39/iK. In the temperature region 
Tc < T < 3Tc, our calculations show that u}zTrf{r = 0) <C 1. Thus the dominant contribution in the integral of Eq. H76|) 
arises from the low density tail of the cloud where ujzTrj{r) ^ 1. In Fig. ^ we plot the temperature dependence of the 
damping rate. For comparison, we also plot the damping calculated using the moment method. We find that the two 
methods give results of the same order of magnitude, but there are significant differences. 



VII. MONOPOLE-QUADRUPOLE MODE IN A BOSE-CONDENSED GAS 

In this section, we consider the m — collective mode in the superfluid phase below Tbec- For this purpose, we 
first need to calculate various equilibrium quantities, solving Eq. (|14|) and Eq. 1)16(1 self-consistently. Fig. [21 shows 
the temperature dependence of the condensate fraction. The effective Bose condensation temperature is lower than 
the free gas result Tc, due to the mean field. Fig. O shows the associated density profile of the condensate and 
noncondensate components at T = 3fiK ~ 0.68Tc. 

Below Tbec J the gas in general exhibits coupled oscillations of the condensate and noncondensate components. At 
finite temperatures slightly below Tbec j one has a "condensate mode" , in which the condensate component mainly 
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FIG. 1: Temperature dependence of the hydrodynamic damping rate of the m — mode in a degenerate Bose gas above Tbec, 
calculated from Eq. (I75II . We also show the result obtained from the moment method by the broken line [see Fig. |HJa)]. For 
comparison, we also plot the moment result for the Maxwell-Boltzmann (MB) gas. 




FIG. 2: Condensate fraction versus temperature. 



oscillates, and a "noncondensate mode", in which the noncondensate component mainly oscillates As one might 
expect, the frequencies of the modes are close to those of a pure condensate mode at T = and a pure noncondensate 
mode above TbeCj respectively. These frequencies are slightly shifted due to coupling between the two components. 
In the calculations in this Section, we focus entirely on the damping of the modes, and neglect these relatively small 
frequency shifts. Thus, we neglect the condensate oscillation in the noncondensate mode and noncondensate oscillation 
in the condensate mode. We calculate the hydrodynamic damping of these two modes. In contrast to the moment 
method results discussed in Appendix B, our present results are only valid in the hydrodynamic regime. 
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FIG. 3: Density profile of tiie condensate and noncondensate along the z axis at T = (~ O.GSTc). The length unit is the 
average harmonic oscillator length ano = h/mu, and the density unit is a^Q- The discontinuous change in rtco and no at the 
condensate boundary is a well-known artifact arising from using the Thomas-Fermi approximation. 



A. Noncondensate mode 



We first consider damping of the noncondensate mode of monopole-quadrupole symmetry. In Eq. H55|) . we use u„ 
given by Eq. (|31() Eq. (|70|l with Eq. H71() and take Uc = 0. For simphcity, we assume that VT = also holds below 
Tbec- We then find that the damping consists of two contributions, F = Fi + F2, where Fi is the contribution from 
the shear viscosity, again given by Eq. (|75(l . while F2 is due to the second viscosity, which is given by 

r. = ^ ^ + (78) 

As noted above, the frequency Q.- of this mode is well approximated by il_ as given in Eq. (|69|l for T > T^- In Fig. 01 
we plot the temperature dependence of the damping of the noncondensate mode below Tbec ■ The contribution from 
the shear viscosity is dominant for T > 0.3Tc, with the contribution from the second viscosity coefficient only taking 
over at very low temperatures. 



B. Condensate mode 



We next consider the condensate mode below the superfluid transition temperature. The pure condensate mode 
frequencies at T = were first given by Stringari [gJl: 



= 2ul + ± i ^16^1-H9^f-16c.>2. (79) 

The associated condensate velocity field is given by 

Uc — {ax,ay,bz), (80) 

where one finds 

a± fnl 3\ h± fnl \ , ^ 



Using this pure condensate mode solution for u^, and setting the noncondensate velocity u„ = and ignoring any 
temperature fluctuation bT ^ we obtain a simple formula from Eqs. Ht)3|) and H34|) for the damping of the low-frequency 
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FIG. 4: Damping rate of the noncondensate mode below the BEG transition temperature. The broken hnes gives the separate 
contributions from the shear viscosity (Fi) and the second viscosity (F2). Gompare with results above Tc shown in Fig. 1. 



condensate mode, 



1 + [yi-T^,Y 

2rn / drricou^ 



(82) 



Here fl- is the frequency given by Eq. H79|l and is approximately fi- « y ^^z, while Uc is given by Eq. H8()(l . both of 
these describing undamped hydrodynamic mode. 

The expression for the damping rate in Eq. (|82|l which has been found here in the hydrodynamic regime is similar 
to Eq. (79) of Ref. l20l|, which gives the coUisional damping of condensate collective modes in the collisionless regime. 
In fact:, one can show that in the limit il-r^ ^ 1, Eq. (|82|l formally reduces precisely to the expression in Eq. (79) of 
Ref. As shown in Ref. this latter ex pres sion for the condensate mode damping is equivalent to the result 

derived by the method of Williams and Griffin '23| in the collisionless regime, under the assumption that the thermal 
cloud always remained in static thermal equilibrium. This makes sense since our assumptions of u„ = and VT = 
are equivalent to assuming a static thermal cloud. 

In Fig. |S1 we plot the temperature dependence of the damping of this condensate mode. The damping of this 
condensate mode is extremely small, simply because we are always in the extreme hydrodynamic limit lUzT/^ <C 1. 



VIII. DISCUSSION AND CONCLUDING REMARKS 

Recently, the Landau-Khalatnikov two-fluid equations were derived ^ trapped Bose gas in the collision- 

dominated local equilibrium domain. These involve the transport coefficients (thermal conductivity, shear viscosity, 
and four second viscosity coefficients) describing the processes leading to local equilibrium in a superfluid, from which 
various relaxation times can be extracted. In the present paper, we used these equations to derive a general 
expression for the damping of hydrodynamic modes in a trapped Bose gas. This formula makes use of the variational 
solution of the two- fluid hydrodynamic equations in the Landau limit {toTi <C 1, where w is the frequency of a 
hydrodynamic mode). We hope that our work will stimulate further experimental studies of the collision-dominated 
hydrodynamic regime in trapped Bose gases. 

As illustration, we used our formalism to evaluate the hydrodynamic damping of the m = monopole-quadrupole 
mode of a cigar-shaped trap, using parameters appropriate to the recent ENS experiments [ifllllllS metastable 
^He*. We presented results for both T > Tbec and T < Tbec- We also used the moment method developed by 
Guery-Odelin et al. j^] for a trapped classical gas and give results (see Appendix B) for a degenerate trapped Bose 
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FIG. 5: Damping rate of the m = condensate mode. 



gas above Tbec- The advantage of the moment approach 19] is that the resulting equations of motion for various 
moments can be solved over the entire frequency domain, including both the coUisionless {lot > 1) and hydrodynamic 
regions {ujt < 1). 

For temperatures characteristic of the ENS experiments (T ~ STc), the effect of Bose statistics is almost negligible. 
Thus, as expected, our calculated values of the damping of the coupled m = monopole-quadrupole mode (due 
entirely to the shear viscosity) are in good agreement with the moment calculations of Ref. jl^] for a classical gas. 
Indeed, our calculation shows that the difference remains small down to the superfluid transition Tc. However, as 
discussed in Ref. |llL , the analysis of the ENS experimental data exhibits a puzzling discrepancy. The damping 
and frequency of the m = mode is consistent with a maximum value of the collision rate [defined in Eq. l)B9|l ] being 
given by Fcou ~ 2 x lO^^s^^, which is achieved at T = 3Tc. Lowering the temperature led to an apparent decrease in 
the values of Fcoii, the latter being determined by the measured changes in the frequency and damping of the ni — Q 
mode. This result seems inconsistent with the calculated value of Fcoii = ^/2nav at Tc, using the measured values 
at Tc of the average density and velocity of atoms, and the collision cross-section a — Sttc^. Our calculated value is 
rcoii(7c) = 9 X lO'^s"^ (see Fig. |7|l, considerably larger than the value Fcou = lO'^s"^ as estimated from the frequency 
and damping of the m = mode [ij [ill • 

Thus, the ENS experiment on the m = collective mode do not seem to be able to enter deeply into the hydro- 
dynamic regime. Ref. ^| tentatively interprets this to being due to increasing inelastic collision processes, which 
effectively lead to a decreasing collision rate FcoU for T below 3Tc. However this does not explain why the estimated 
value of FcoU calculated at Tc is an order of magnitude larger, which would correspond to the m = mode being 
deeply into the collision-dominated hydrodynamic domain. 

In order to clarify this puzzling behavior above Tc, it would be useful to measure the damping and frequency of the 
monopole-quadrupole collective mode in the Bose-condensed region. As discussed in Ref. the effective value of 
the relaxation time associated with the shear viscosity is calculated to become much smaller as one goes below Tc. 
This is simply because l/r^ of the thermal cloud atoms is dominated by coUisions with condensate atoms. As a result, 
even if one has a low thermal cloud density (spatially averaged) such that one is in the coUisionless region above Tc, 
one is automatically deep inside the hydrodynamic region for T < Tc because of the rapid build up of the condensate 
density in the center of the trap. This suggested variant of the ENS experiments effectively uses the formation of 
the high-density Bose condensate to increase the collision rate which determines the frequency and damping of the 
monopole-quadrupole mode. 

In Section VII, we presented the first explicit calculations of the hydrodynamic damping of the monopole-quadrupole 
mode in the superfluid phase. As can be seen in Fig.^ the damping of this mode involving the noncondensate thermal 
cloud is still dominated by the shear viscosity down to T ^ 0.3Tc. Comparing the results above (Fig. ^ and below 
Tc (Fig.^, one sees the mode damping F is fairly smooth going through the transition, with only a slight decrease in 
magnitude below Tc- This slight decrease in the value of F hides the fact that (as discussed above) the effective value of 
the shear collision relaxation time r is rapidly decreasing as we go below Tc, putting one deep into the hydrodynamic 
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domain. 

In Section VII, we also evaluated the hydrodynamic damping of the monopole-quadrupole mode in the condensate. 
As seen in Fig. [3 these damping is extremely small since one is effectively in the "Landau limit", lot^ <C 1. 
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APPENDIX A: FREQUENCY-DEPENDENT TRANSPORT COEFFICIENTS 

In this Appendix, we give some details on the frequency-dependent coefficients of the shear viscosity and thermal 
conductivity, starting from the kinetic equation for the noncondensate atoms. As in the usual Chapman-Enskog 
procedure described in Rcf . P| , we insert the local equilibrium distribution function /igq in the left hand side of the 
kinetic equation, where /icq is given by 

/lcq(l", P,i) = ^-l^j. ^■jg/3(r,t)[p-mv„(r,t)]V2in _ 2' ^^^^ 

where z{y,t) = e''^'"'*) is the local fugacity. As shown in Appendix A of Ref. 'E^, the kinetic equation is 
then given by 



5/ 

dt 



_i_ . Vz + — . vr + • ( — • V 

z m 2mkBT'^ ni ksT \m 



— - • (p-TOV„) 



/lcq(l+/lcq) =Cl2+C22. (A2) 



In contrast to Ref. ^ , we keep the time derivative of / explicitly. Since we are interested in small-amplitude collective 
oscillations, in the following we always expand the theory to first order in the fluctuations around static equilibrium . 

We first consider the shear viscosity, which is associated with the anisotropic pressure tensor. In a linearized theory, 
this is given by 

PI, = P^, - S^,P = J (p^P. - y,-P') /(r, P, t)- (A3) 

The equation of motion for P^' ^ can be obtained by taking moment of Eq. (|A2|) and linearizing it around static thermal 
equilibrium. One finds that the term contributes to this moment is 

dp 1 f 1 r 2\ P fP ^\ f n , f \ - 5 f ^'""'^ , ^""P 2 



One thus obtains 



where 



^or^ + ^-^5..V.v„^ - . (A5) 



dt \ Br Br Q-t^" ' \r7i 3 



coll 



coll 



— PiJ^Pu - TlS^uP ]) = / 77r-rvT~ Pi'P'' ~ o^l^uP (C12 -I- C22). (A6) 



The collisional contribution on the right hand side of Eq. IIASII arises from deviation of the distribution / from the 
local equilibrium solution in Eq. (|A1() . Following the Chapman-Enskog procedure, we use the ansatz / = /icq + (5/, 
where 

5f = Y, B^, fp^p, - ls,,pA /o(l + /o), (A7) 
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with Bfj^iy being some momentum-independent symmetric tensor. The relation between B^i, and P'^^ can be found by 
using Eq. (|A7() in Eq. (|A3() and carrying out the momentum integrah 



dp 1 



1 



3 

dp 1 



1 



h)^ m 



Pt^'Pi'' - o Vi^'-P' 



Using Eq. (HT)) in Eq. (|X6|) . we find 



- Up. - 3<5,.p- 



dp 1 



(27r;i)3 TO 



1 



(A8) 



(A9) 



where L is the hnearized coUision operator defined in Eqs. (46) and (48) of Ref. Combining Eq. (|A8(I and Eq. (|A9p 
and using the definition of r,, given in Eq. (B5) of Ref. |^], we find that the coUision term reduces to 



- Up. - ^s,.P 



p' 



(AlO) 



where the viscous relaxation time t,, is defined in Refs. j^l^' Assuming the harmonic time dependence P^j, cx e 
we finally obtain 



P' 



2t,,Po 

1 — iLUT„ 



(All) 



where the frequency-dependent viscosity coefficient 77(ti-') is defined by 



T]{UJ) = 



1 — ilUTrf 1 — ilOTrf 



(A12) 



The frequency-dependent thermal conductivity can also be obtain in the same manner, by considering the linearized 
heat current 



Q(r,i) = 



P^ V> rr .^ 5 ~ 



dp 



p-" 5 g5/2(^o) 

2m 2 53/2(20) 



-/(r,P,i). 

m 



Taking the moment of the kinetic equation in Eq. |IA2|I . one finds that the relevant contribution is given by 



(A13) 



dp 



{2TThf 
2 TO 



5 55/2(^0) 
2m 2 93/2(^0 j 



V6T 



/ 757/2(^0) 
\ 233/2(20) 



TO 2mk^TQ Vto 
55/2(20)"'^ 



(P .V<5r)/o(l + /o) 
Vto / 



53/2(20) 



One thus obtains 

aq 5 fiQklTQ 



dt 



VST 



/ 757/2(20) 
\ 253/2(20) 



where 



5 .g5/2(2o) 

feBJo — ^ — 7 T 

2m 2 53/2(20) 



5 


55/2(20)" 


X 


2 


53/2(20), 








" 5 


~ J 


' (27r?i)3 


2to ^ 2 



5 55/2(20) 

kbJo — — } — 7 

2to 2 53/2(20) 



«B-'0 



55/2l2oj 
53/2(20) 



— (C12 + C22) 
TO 



(A14) 



(A15) 



(A16) 
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To evaluate the collisional term Eq. (IA16|I . which arises from deviation from local equilibrium, we use the Chapman- 
Enskog ansatz / — /icq + Sf, where j5|| 



5/ = A. 



5 55/2(^0) 

fceJo — — 7 — r 

2™ 2 53/2(2:0) 



/o(l + /o). 



Here A is a momentum-independent vector which is directly related to the heat current Q through 

\-i 2 



Q = A 



dp \ p^ 



{2TTnY 3m 



5 ^5/2(^0) 

7, oKbJo 7 T 

2m 2 53/2(2^0) 



/o(l + /o). 



(A17) 



(A18) 



Evaluation of the collisional term with using the ansatz Eq. (|A17|) closely follows the derivation of the thermal 
conductivity in Ref. (5j . We find 



rp2 



5 55/2(^0) 
2m 2 53/2(^0) 
2 5 



A 

y 



r t), ^ 55/2(2:0) 

-/cbTo — — 

2m 2 53/2(20) 



L 



P^ 5 55/2(2:0) \ p 

Kb^O — — 7 T (■ — 

2m 2 53/2(20) J m 



(A19) 



Combining Eq. (|A18|I and Eq. (|A19p and using the definition of the thermal relaxation time given in Eq. (Bl) of 
Ref. 0, we find 



5 55/2(20) 
2m 2 53/2(2:0) 



Assuming the harmonic time dependence Q cx e we obtain 



coll 



9l 



Q 



1 - ZWTk 



V6T = -k{uj)VST, 



(A20) 



(A21) 



where we have used the expression for k given in Eq. (|24|) . and the frequency-dependent thermal conductivity is 
defined by 



1 - icjTK-, 



(A22) 



APPENDIX B: MOMENT METHOD FOR A DEGENERATE NORMAL BOSE GAS 

In the moment method |l9j |. one derives an equation of motion for some dynamical quantity denoted by x(r,t) by 
taking a moment of the kinetic equation in both position and momentum: 

drj^x{r,mr,p,t). (Bl) 

The advantage of this method is that it gives results valid in both the coUisionless and hydrodynamic region. It should 
be noted that in the coUisionless region, this approach does not include Landau damping, but this is small above 
Tbec- For the m = monopole-quadrupole mode, we need moment equations for the following physical quantities: 

Xi = X2 = 2z2-r2, 

X3 = r-p/m, X4 = 2zpz/m - r_L • p^/m, 

X5=pVm^ X6^2pl/m^-pl/m\ (B2) 
Calculation of these moments give the following coupled equations 

^-2(X3) = 0, fel-2(x4)^0, 

^-(X5) + ^^^(Xl) + ^^(X2)=0, 
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FIG. 6: Moment calculation of the (a) frequency and (b) damping of the m = mode, both as a function of the quadrupole 
relaxation time r. These results are found by solving Eq. I|B11^ . We use the ENS trap frequencies ljx/Stt — 988Hz and 
Lj,/27r = 115Hz [l3[ll[l3- 



(X3 



0, 



dt 3 

— ^ ^ U3/ H ^ \X4/ — U6/coU- 

The collisional contribution in the last equation in ljB3(l is defined as 

This term can be approximately evaluated by using the ansatz f ~ fo + Sf, where 
'5/-/3o/o(l + /o) 



ST ( 

-t^\^^Uq- ficQ ) + p • v„ + + a{2pl - 

1q \ Z7TI 



(B3) 



(B4) 



(B5) 



where the time-dependent parameter a{t) characterizes the anisotropy in the momentum distribution described by 
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FIG. 7: Temperature dependence of the spatially-averaged quadrupole relaxation time r. The broken line shows the result 
using the Maxwell-Boltzmann distribution for a classical trapped gas. 



(Xe)- The above ansatz is a generalization of the gaussian ansatz for the Maxwell-Boltzmann gas used in Ref. [T^ to 
a degenerate Bose gas. 

Using the ansatz in Eq. (|B5|I in Eq. HB4|I and linearizing in a, one obtains 

(X6)con = (B6) 

T 

where r is a quadrupole relaxation time defined by the weighted spatial average of the inverse of the viscous relaxation 
time (see also Ref. l2(|) 

-r JdrPa 



Using the moment equations Eq. IjBSp together with the approximation Eq. IIB6|) , the only effect of using Bose statistics 
is in the value for the spatially-averaged relaxation time r in Eq. 1B7() . 

In a nondegenerate (Maxwell-Boltzmann) gas, one finds that r^(r) — |Tci(r) [ElEll^i where t^i is the usual elastic 
collision time for a classical gas 

r^7\r) = \/2o-no(r)5 = %/2(87ra2)no(r)(8fcBT/7rm)i/2_ (gg) 

Using this result and Po(r) = kBTno{r) in Eq. HB7|I , the quadrupole relaxation time reduces to 

1 4 

- = ^rcou. (B9) 
r 5 

Here Fcoii = \/2afiv, where the spatially-averaged density is defined by 

no(r = 0) 



drng(r) 



(BIO) 



This result for l/r in terms of Fcoii agrees with that of Guery-Odelin et al. . The authors of Ref. fT3| analyzed 
data in terms of the classical gas result of Ref. |l^J , using the above definition of the averaged collision rate Fcou valid 
for a Maxwell-Boltzmann gas. 
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Assuming the time dependence e , the coupled equations in Eq. IjBSp with Eq. ljB6|l can be solved to give 



0. 



(Bll) 



Solving Eq. (|B11|) . we obtain the solution u)^Vl — iT, describing damped modes. One can see that in the collisionless 
limit UJT ^ 1, the frequency is given by either 2ljJz or 2a; j^. In the opposite hydrodynamic limit lot <^ 1, the two 
solutions are given by w = fl±, with fl± given by Eq. H69|l . The frequency fl and damping T obtained by solving 
Eq. (|B11I) can be expressed in terms of the spatially-averaged quadrupole relaxation time r, for given trap frequencies. 
In Fig. El we plot fl and T of the low-frequency mode as a function of ujzT. The hydrodynamic domain is the region 

LOzT < 1. 

In Fig. [7| we plot the temperature dependence of the quadrupole relaxation time t calculated for the ENS experi- 
mental data. For comparison, we also plot the result [i3,[l3using the Maxwell-Boltzmann distribution. The effect of 
Bose statistics is clearly very small, down to about T ~ 1.5Tc. We use the result in Fig. [Tito calculate the frequency 
and damping as functions of the temperature, as shown in Fig. 

In the hydrodynamic limit, one can obtain an analytical expression for the damping rate from the moment equations. 
Assuming ojt ^ 1, one can expand the solution to Eq. ljBll|l to first order in r. In this limit, the damping r_ of the 
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low-frequency mode fi^ is given by 

r- = - ^^im'- - (B12) 

Apart from a different averaged shear-viscous relaxation time r, it is satisfying that this moment result is identical 
to the LL expression in Ea. lf75)l . 
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